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Abstract. We study the impact of Early Dark Energy fluctuations in the linear and non- 
linear regimes of structure formation. In these models the energy density of dark energy is 
non-negligible at high redshifts and the fluctuations in the dark energy component can have 
the same order of magnitude of dark matter fluctuations. Since two basic approximations 
usually taken in the standard scenario of quintessence models, that both dark energy den- 
sity during the matter dominated period and dark energy fluctuations on small scales are 
negligible, are not valid in such models, we first study approximate analytical solutions for 
dark matter and dark energy perturbations in the linear regime. This study is helpful to 
find consistent initial conditions for the system of equations and to analytically understand 
the effects of Early Dark Energy and its fiuctuations, which are also verified numerically. In 
the linear regime we compute the matter growth and variation of the gravitational poten- 
tial associated with the Integrated Sachs- Wolf effect, showing that these observables present 
important modifications due to Early Dark Energy fluctuations, though making them more 
similar to ACDM model. We also make use of the Spherical Collapse model to study the 
influence of Early Dark Energy fluctuations in the nonlinear regime of structure formation, 
especially on 6c parameter, and their contribution to the halo mass, which we show can be 
of the order of 10%. We finally compute how the number density of halos is modified in 
comparison to ACDM model and address the problem of how to correct the mass function 
in order to take into account the contribution of clustered dark energy. We conclude that 
the inhomogeneous Early Dark Energy models are more similar to ACDM model than its 
homogeneous counterparts. 
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1 Introduction 

The understanding of the accelerated expansion of the universe is one of the greatest chal- 
lenges in physics. If we assume that General Relativity is our accepted theory of gravitation 
we must introduce a new form of fluid with sufficiently negative pressure, p < — /o/3, that 
accounts for roughly 3/4 of the universe energy density today. The physical description of 
this new form of matter, generically called dark energy (DE), is yet unknown. Latest data 
analysis of luminosity distance of Supernovae type la [1], cosmic microwave background [2, 3] 
and large scale structure [4, 5] are all consistent with a flat universe with approximately 1 /4 
of its critical density in the form of pressureless matter (cold dark matter and baryons) and 
3/4 in the Cosmological Constant, A. 

If the accelerated expansion is caused by A, since it is constant in space and time, 
it does not cluster and has a negligible contribution to the energy density budget of the 
universe at high redshifts, affecting solely the background evolution for z ~ 1 and lower. 
Although A is the simplest model for the accelerated expansion, it suffers from two severe 
theoretical problems. Since the natural interpretation of A is the vacuum energy, Quantum 
Field Theory should be able to determine its value. However it predicts a value for A that 
can be several tens of orders of magnitude larger than what is observed, which is known as 
the Cosmological Constant Problem, see, e.g., [6, 7]. Another problem, more closely related 
to the cosmological evolution itself, is why the observed value of A is such that it becomes 
important for the evolution of the universe just at the time we are able to measure its effects, 
which is known as Coincidence Problem [8]. 

Many alternative models for the accelerated expansion have been proposed. Possibly 
the most studied ones are based on canonical scalar fields [9-11], which are usually referred as 
quintessence models. The evolution of quintessence is background dependent, fact that could 
explain the transition from a decelerated to an accelerated phase as a natural evolution 
between attractor regimes [8, 12]. This kind of mechanism can potentially alleviate the 
Coincidence Problem, diminish the dependence on the initial conditions of the scalar field 
and bring new features to the cosmological evolution, for instance, the possibility that the 
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energy density of DE is non-negligible at high redshifts. Models presenting this behaviour 
are called Early Dark Energy (EDE) and were extensively studied in the literature, see for 
instance [13-17]. 

Another important feature of dynamical DE models is that, in contrast with A, they 
possess fluctuations. On small scales, in the linear regime, quintessence perturbations are 
several orders of magnitude smaller than Dark Matter (DM) perturbations and are usually 
neglected in studies of structure formation. This is due to the fact that the effective sound 
speed of canonical scalar fields perturbations, or the sound speed in the rest frame [18], 
is CefT = 6pc/6pc = 1, which suppresses the growth of field perturbations inside the sound 
horizon scale, which, in turn, is of order of the particle horizon. However quintessence 
fluctuations can not be neglected from both the theoretical and observational point of view 
[19, 20]. 

Moreover, there exists some realisations of DE models where its fluctuations can grow on 
sub-horizon scales, i.e., with Cggf <C 1, such as k-essence models [21-24]. The possibility that 
DE has an effective sound speed less than unity has been investigated by many authors, e.g., 
[25-29]. In particular, Ref. [29] points out that current CMB and LLS data slightly prefers 
dynamical DE, Ces ^ 1 and some amount of EDE. It is also worth to note that clustered DE 
seems to give a better prediction for the concentration parameter of massive galaxy clusters 
[30]. 

Structure formation in EDE models has been studied by many authors, e.g., [14, 15, 31- 
36]. In particular, for the case of Ces = 1, Ref. [34] shows that neglecting EDE perturbations 
leads to incorrect constrains on the equation of state and Ref. [35] claims that EDE models 
can be constrained by future observations of galaxy clusters. However all these studies either 
neglected EDE perturbations or consider models with c^s = 1, which effectively renders 
negligible perturbations on small scales. 

The objective of this paper is to analyse structure formation in EDE models that can 
present large fluctuations on small scales. In this scenario DE has two major characteristics 
not present at the same time in the usual quintessence models: 1) DE energy density is 
non-negligible at high— z and 2) DE fluctuations can be of the same order of magnitude of 
DM fluctuations. We compare the results with the usual assumption of nearly homogeneous 
EDE and ACDM model. Assuming that EDE is described by a perfect fluid, characterised 
by its equations of state, pe = w{t)pe, and the effective sound speed of its perturbations 
Cgg = 6pe/5pe, we analyse both the linear and nonlinear evolution of EDE fluctuations and 
their impact on DM growth, compute the number density of halos and how it is modified by 
the contribution of DE fluctuations. 

The outline of the paper is the following. In Sect. 2 we present and discuss the back- 
ground evolution of two models of EDE. In Sect. 3 we study the evolution of linear pertur- 
bations of a system with EDE and pressureless matter and calculate the matter growth and 
the Integrated Sachs- Wolf effect. Sect. 4 is devoted to the study of the nonlinear evolution 
and the Spherical Collapse Model. In Sect. 5 we present the results for mass functions and 
in Sect. 6 we present our conclusions. 
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2 Background evolution 



We assume a universe with flat spatial section, DM (baryons are treated as dark matter) and 
DE. Friedmann's equations in conformal time are then: 

= ^a^^Pm + Pe) and n = -^a^[pra + Peil + 3w)] , (2.1) 

where Ti = a/a and the dots represent derivative with respect to conformal time. We choose 
two different parametrizations for EDE. In the first one the energy density parameter of DE, 
Qe = 87rGa^/?e/3?^^, is given as a function of its value at early times, its equation-of-state 
parameter now, wq, and matter energy density now, [13]: 

K «l(l-a-^-). (2.2) 

We will refer to this parametrization as Model A. The advantage of this parametrization is 
that one can fix fi^, however it does not allow to directly choose details of evolution of the 
equation-of-state parameter. In order to control some properties of w, such as the value in 
the matter dominated era, w^a, the moment of transition from w^a to wq, Cc, and the duration 
of this transition, Am, we also study the following parametrization [37]: 

w (a) =wo + [Wm - Wo) -. ^ r-TV ' (2-^) 

We refer to this parametrization as Model B. Using the parametrization (2.3), on the other 
hand, we need to fix values of the parameters Am, Oc and in order to give the desired 
value of Og. In this case the DE energy density is given by: 

n f r {l + wia'))da'\ , 
p, (a) = pI exp (^-3 ^ ^ J . (2.4) 

For both models we choose the amount of DE at early times to be ^% ~ 0.018, consistent 
with the limits presented in [17], the amounts of matter and DE now are f^m — 0-25 and 

= 0.75 and the DE equation of state now is = —0.9. In Model B we also set w^a = 
—0.1655, ac = 0.5 and Am = 0.09. For these two models we show the evolution of f^m (o), 

(a), w (a) in Fig. 1. 

As we can see in Fig. 1, in Model A, ilg (o) is basically constant at high-z and its 
equation of state varies slowly, whereas in Model B the amount of DE has a non-negligible 
variation during most of the cosmic time and has a rapid transition of its equation of state 
for o ~ 1. We will show that these differences in the background evolution will imprint 
distinct features both in matter and in DE fluctuations. Assuming a Hubble constant of 
= 72kms"^Mpc"^ the age of the universe in a ACDM model with = 0.25 is 13.77 Gy, 
whereas in Model A we have 13.44 Gy and in Model B 13.01 Gy. 

In Fig. 2 we also show the evolution of the comoving volume, given by 



(2.5) 



dzd^ H {z) ' 

where r {z) = {z')dz'. This quantity, which depends only on the background evo- 

lution is important for the study of cluster number counts, which in turn also depends on 
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Model A 



Model B 




a" 



Figure 1. Evolution of matter and DE density parameters, J7m (a) (solid black line) and fide (a) (blue 
dashed line), DE equation of state w (a) (red dot-dashed line). Left (right) panel: Model A (B). 
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Figure 2. Comoving volume for ACDM model (solid black line) and models A (blue dotted line) and 
B (red dashed hue) of EDE. 

perturbative properties via the mass function. Note that both models of EDE that we are 
considering present a smaller volume then ACDM model and that important differences ap- 
pear only at high— z. For redshifts z < 0.5 all three volumes are very similar, which already 
suggests that cluster observations at low redshifts would poorly differentiate between these 
DE models. We will return to this issue in Sect. 5. 



3 Linear evolution 



In this section we study the linear relativistic evolution of a system with DE and matter. In 
models of EDE the energy density of DE at high redshifts, e.g., at the redshift of decoupling 
Zdec) is not negligible, hence we carefully analyse approximate analytical solutions in order to 
establish consistent initial conditions for the equations of motion. Moreover this study will 
clarify some effects due to EDE and its perturbations. 

In the Newtonian gauge, in Fourier space, in the absence of anisotropic stress, Ref. [38] , 
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the perturbed equations as a function of conformal time can be written as: 



where 



^m + ^m = 3(/., (3.1) 
ern+nem = k^<P , (3.2) 

(5e + m (^j^ -w^5e + {l + w)e^ = Z{l + w)^, (3.3) 

e. + n{i- 3cl) = ^^^f^^^ + k'cp , (3.4) 

[1 + w) 

k^cj) + 3n[^ + n(j)) = — — {n^6,^ + rieSe) , (3.5) 

2 Pe W 
Ca = — = w — — r (3.6) 

is the squared adiabatic sound speed and the pressure perturbation is given by [26]: 

^Pe = clffdpe + 3n{l + W) [els - Ca) Pe^ ■ (3.7) 

Note that the equation of state can not cross the phantom barrier, otherwise wiU diverge, 
and so will DE perturbations. In the two models we study there is no phantom crossing. 
For a treatment of this issue see Ref. [39]. We solve the system of equations (3.1)-(3.5) 
numerically, however the (ii)-Einstein equation, given by 

'<p + 3n^+(--n'\^= (p)se, (3.8) 



2 \Spe 

will be useful to find analytical solutions. 

It will also be useful to write a second order equation for the density contrast. For the 
sake of simplicity we assume that Ceg and w are constants, then we have: 

6 + [n{l-3w) - A]S + 3A (n + - S-HAA + clgk^ 6 = {l + w){S- k^<p) , (3.9) 
where A = c^g — w, 

^-m^A + k^ ^^-^"^ 

and 

S = 34> + 3'H{l- 3cls) 4> - 3A^ - m'^Acf) . (3.11) 

Eq. (3.9) is valid for any perfect fluid characterised by constant w and c^g, so the correspond- 
ing equation for matter is given by assuming w = c^g = 0. 

The function A is related to the presence of the pressure perturbations that are pro- 
portional to 0. It is useful to note that on small scales. A: » 7^, >1 ~ and 6pe — c^^Spe, 
as already observed in [40]. On large scales, k <^7i, we have A ~ 2'H/'H and the pressure 
perturbation strongly deviates from the value prescribed by Cgff . 

There are two important scales that determine the qualitatively behaviour of DE per- 
turbations: the particle horizon 

da' 



>^H{a)= I — , (3.12) 



- 5 - 



(horizon for short), which is also important for matter perturbations, and the sound horizon 



. f s f Ccs (a) da 

As (a) = / • (3.13) 

Perturbations with wavelength smaller than Xg oscillate with decreasing amplitude and even- 
tually reach a minimum value proportional to the gravitational potential. Perturbations with 
wavelength larger than As effectively behave as pressureless, growing at the same pace as mat- 
ter perturbations. Finally perturbations in both matter and DE with wavelength larger than 
Xfj will follow the time dependence of the gravitational potential, being constant whenever 
DE effect in the background is negligible. 

We are interested in small scales were nonlinear structures, such as galaxy clusters, 
form. Typically we can assume that the order of magnitude of such scales is 10 h~^ Mpc, or 
/cni ~ 0.63hMpc^^. Shortly after decoupling, e.g., at z = 1000, the horizon wave number 
is kn — 0.017hMpc"^ (ACDM with QjJi = 0-25 background is assumed), hence both DE 
and DM perturbations that go nonlinear are well inside the horizon on the onset of matter 
dominated period. For this reason we will focus our analysis on scales much smaller than 
\h, consequently the only scale important for DE perturbations is A^. For a study of scales 
larger than the horizon see Ref. [40]. 

It is interesting to determine the value of Ces, assumed constant, such that kg = k^y 
today: 

c„i = 9.4 X 10"^ . (3.14) 

In models with Cgfr <C Cni the sound horizon is always smaller than the nonlinear scale, then 
DE perturbations will behave as pressureless and we can effectively assume Ces = 0. On 
the other hand, in models with Ces 3> Cni, DE perturbations are suppressed by its pressure 
support on nonlinear scales. In this work we treat the two limiting cases: one with negligible 
sound speed, Ces = 0, and one with a non- negligible sound speed, Ces = 1. For studies of time 
dependent Cgff in linear theory see Ref. [41] and for arbitrary values of Ceg during nonlinear 
evolution of DM see Ref. [42, 43]. 

3.1 Non- negligible Ccff 

On scales well inside the horizon, for constant Ces and w, from Eq. (3.9), DE perturbations 
obey the following equation: 

^ ^ a V 'H^ J 
where ' = d/da, Wt = ^eW, 

3 

a = - (1 — Wt) — Sw , and /3 = 
Eq. (3.15) has the particular approximate solution: 

. ^, (1 + ^)0 



which becomes more accurate the bigger Cesk/Ti is. Note that since Cegk/'H grows in time, 
as long as the background parameters do not change too fast in time, <5e will approach 



(1 + w) k'^(j) 



3A 



1 - 3wt) 



(3.15) 



(3.16) 



(3.17) 
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solution (3.17) on all scales. Moreover, the homogeneous solutions of Eq. (3.15) oscillate 
with decreasing amplitude [40, 44], hence eventually solution (3.17) will be reached and can 
be used to set the initial value of 5e as a function of (j). 

Before we determine the initial value of ^e, we proceed to find solutions for (j). Since we 
are interested on small scales we can assume that 6pe — Cg^Jpe- We verified that, for the IC 
we will determine and c^q = 1, the error in this expression compared to Eq. (3.7) for a scale 
k = O.lhMpc"^ is initially of 1 - 3% and below 1% for z < 100 and that the error decreases 
for smaller scales. Then, substituting the expression (3.17) in Eq. (3.8) and changing the 
independent variable to the scale factor a, the evolution of the gravitational potential is given 
by the equation: 

This equation is valid as long as 0,^ and w are nearly constant, which is true for Model A 
roughly for 10 < z < 1000 and for Model B roughly for 100 < z < 1000. Solving this equation 
in this interval enables us to analytically evaluate the time variation of (p. The solutions of 
Eq. (3.18) are power laws (p on a", where n is solution of the algebraic equation 

2 /5 3 \ 3{l + w)ne 



n'+\^---wtjn+^ ^-^^-3u;t = 0. (3.19) 

For non-EDE models, when DE is very subdominant, we can set ilg = and the solutions are 
n = and n = —2.5, which indicates that, neglecting the decaying mode, the gravitational 
potential is constant during matter dominated era, a well known result for the Einstein-de- 
Sitter universe. However, in EDE models (p is not constant any more, its early time variation 
depends on both ile and w, but note it is independent of Ces- 

The largest values of n, which give the solutions that decay more slowly, are 

rimod A = -0.0100 and n^od B = -0.0127. (3.20) 

Since in Model A both Qc and w are approximately constant before DE domination, the error 
in this expression against the full numerical solution is below 1% in the range 5 < z < 1000. 
Given that in Model B the time variation of ilc and w is larger, its analytical solution has 
an error below 1% only in the range 150 < z < 1000. 

Once we know how (p initially varies in time we can find 6c from Eq. (3.17) and use 
Eq. (3.3) to determine the initial value of Be- 

,^^^_l_m^^, (3.21) 

Hence Eq. (3.21) can be used to determine the initial value of as a function of the initial 
value of (p and its initial time derivative. 

Finally we look for analytical solutions for matter perturbations in order to determine 
IC for (5in and Oya- Setting w = c^s = in Eq. (3.15) we have: 
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It is useful to note that, although in EDE models Oe is non-negligible at high-z, on small 
scales, since <5e ~ we can assume that ilm'^m 3> r^e^e in Eq. (3.5), then we have: 

" 3^ • ^'-''^ 

Using Eq. (3.23) we can find the initial value of 6^1 as a function of the initial value of a given 
(j). Inserting Eq. (3.23) in Eq. (3.22) we have the equation for 6^ alone: 



/ 3 3u;t\ 30m 5m „ 



which is the usual equation for the evolution of matter perturbations well inside the horizon. 
The solutions are power laws oc a^, where p is determined by: 

Note that, under the approximations taken, when we compare Eqs. (3.25) and (3.19), we find 
that n = p — 1. Then the corresponding values for p are: 

Pmod A = 0.9900 and p^od B = 0.9873 . (3.26) 

In order to find the initial value of 6^ we can use Eq. (3.1). The initial value of 6m and 
its initial dependence with a, is given by finding p in Eq. (3.25). Then on small scales we 
have: 

Bra = -Sm ■ (3.27) 

In Fig. 3 we show the evolution of 6c according to the numerical solution of equa- 
tions (3.1) — (3.5), and the analytical solution of Eq. (3.17), with </> given by the numerical so- 
lution. The only value we need to choose is (pi = (f> (ai), taking Oj = 0.001 and (pi = —6.2 x 10~^ 
we get 6m (1) — 0.1. In order to check the consistency of the analytical solutions that we 
have used to determine the IC of DE perturbations, given by Eq. (3.17) and Eq. (3.21), (red 
dotted-dashed line) we also show the evolution of 6e using 6e ten times greater than as given 
by Eq. (3.21) (blue line). As we see, the solution of Eq. (3.17) is in good agreement with the 
numerical solution. In the case of bigger 9^, 6e oscillates with an initial greater amplitude, 
which decreases with time and eventually approaches the value given by Eq. (3.17). For 
sake of clarity in these plots, we assumed c^g = 0.1 in order to have a smaller frequency of 
oscillation. It is important to note that for this scale and for this value of Ces, at low z, DE 
perturbations are 4 orders of magnitude smaller then 6m- For smaller scales and larger Cgg, 
6p will be even smaller than 

3.2 Negligible CeS 

As we will see, in the case of negligible Cefj, DE perturbations can have the same order of 
magnitude of matter perturbations. Although this fact impedes us to assume that r^m'^m ^ 
fte6e for EDE models, fact that simplified the equations for the case of non- negligible c^s, 
when Ceff is negligible we will have another sort of simplification. Note that substituting the 
expression for the pressure perturbation, Eq. (3.7), in the equation for 9^, Eq. (3.4), and 
assuming Cgfj = we obtain 
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b.ool 0.01 0.1 I 0.001 0.01 0.1 1 

a a 

Figure 3. Linear evolution of 5^ with c^jj = 0.1 for a mode k = O.lhMpc^^: analytical solution, 
Eq. (3.17) (black dashed line), numerical solution with IC from Eqs. (3.17) and (3.21) (red dot-dashed 
line) and the numerical solution with O^, ten times bigger (blue line). The change in the analytical 
solution due to the different IC is almost imperceptible in this scale so we only show the first case. 



e^+n9^ = k^(t), (3.28) 

which is just the same equation for ^m, Eq. (3.2). Hence both matter and DE perturbations 
will feel the same force and flow in the same way. Then, using Eqs. (3.1) and (3.3), we can 
determine an equation that directly relates and 5e'- 

5e — SHwde = [1 + w) 6ra , 

where, since, on small scales, both 6^^ and 6c will be of order 'H~'^k'^4', we neglected the 
(f> term. Assuming that oc ci^ we find the following relation between DE and matter 
perturbations: 

= ^^5. . (3.29) 

p — 611! 

In models with negligible amount of DE during matter era we have p = 1, then 



(1 - 3w) 



which is in accordance with literature [28, 44, 45]. However, as we will show in Eq. (3.34), 
p ~ 1 in EDE models with negligible c^s, then we can actually use Eq. (3.30) to set the initial 
value of de as a function of dm- It is also interesting to observe that Eq. (3.30) indicates 
that DE perturbations diverge as w — )• 1/3, which can be achieved by a tracking scalar field 
during radiation dominated era. However this solution is valid only for Ces = 0, hence, as 
consistency condition for DE perturbations, we should expect that whenever w ^ 1/3 we 
must have Ces > 0. In the models we are studying, during the matter dominated era and DE 
domination, we never have w = 1/3. 

At this point we derived a relation between DE and matter perturbations but not their 
evolution with time. Setting w = Ces = in Eq. (3.15) we have the equation for the evolution 
of (5m: 



< + ( ? - ^) ^ = ^ i^mSm + njc) . (3.31) 
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Now we can use Eq. (3.29) to express 5e as a function of 6^1 and rewrite Eq. (3.31) as 

While and w are nearly constants we can find the solution tx a^, where p is 
solution of 

,.,,(^_^)_^(a.,.<l±|4.„.).„. (3.33, 

For the two models we are studying we have the following values at af. 

Pmod A = 1.0052 and p^od B = 0.9934 . (3.34) 

As we can see, the deviation from p = 1 is very small and we could actually use expression 
(3.30) in order to determine the initial value of 5e as a function of d^a- Moreover expression 
(3.30) gives a better order of magnitude estimation of 6c for low— 2; than Eq. (3.29), that 
depends on the determination of p from Eq. (3.33), which is not accurate for low— z, when 
both and w can not be assumed constant. Hence in what follows we will use the expression 
(3.30) to compute the initial value of the DE perturbation. 

To determine the initial value of 5^ as a function of (f) we use Eq. (3.29) in Eq. (3.5), 
then, for small scales we have: 

(3.35) 



37^2 (n^ + n^^ 



3w 

The initial value of can be determined in the same way as we did for the case of non- 
negligible Ceff, using Eq. (3.27). 

Now we proceed to find the initial time variation for <j). For the case Ces = DE 
perturbations do not contribute to the variation of (p, i.e., the RHS of Eq. (3.8) is zero and 
then the equation that determines the power of n of (/> is the following: 

2 /5 3 



n' + n\^-- -wtj -3wt = 0. (3.36) 

Here, in contrast to the case of non- negligible Ces, the relation n = p — 1 is not valid. Note 
also that n only depends on the product OgW^, differently from case of non- negligible Cgg, 
where it also depends directly on ilg. This indicates that EDE models with negligible sound 
speed give a very low contribution to the time variation of due to the fact that, although 
Qe is not negligible during the matter era, we have — 0. The highest values of n for 
models A and B are: 

'^mod A = 0.0017 and n^aod B = —0.0036 . (3.37) 

Note that the absolute value of n is one order of magnitude smaller than the solutions for the 
case of non-negligible Ceff, Eq. (3.20). Since the background evolution is the same for both 
cases, this clearly shows that the most important effect for the evolution of the gravitational 
potential in EDE models is pressure perturbation. 

In Fig. 4 we show the evolution of matter and DE perturbations as well as the evolution 
of 6e according to Eq. (3.30). We show this relation, instead of the more general one, in 
Eq. (3.29), because initially the deviation from p = 1 is very small (bellow 1%) and for late 
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Model A 
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Figure 4. Linear evolution of perturbations with Ccff = and k = O.lhMpc^^: numerical solution 
of (5ni (black solid line), numerical solution of (blue dashed line) and dc according to the relation 
given by the analytical solution, Eq. (3.30), (red dot-dashed line). 

times, when Jim, f^e and w are not constant, the computation of p via Eq. (3.33) is less 
accurate than just assuming p = 1. As we can see, Eq. (3.30) is a very good approximation 
for 6e during the matter dominated era, however, once the transition to the DE dominated 
era starts, it becomes much less accurate. Anyhow Eq. (3.30) can be used as an estimate 
of the order of magnitude of 5^ for late times. It is also interesting to note that in Model 
A the perturbations in DE are larger than matter perturbations until w starts to decrease, 
a ~ 0.1. For Model B, since 1 + it; is smaller than in Model A, 6c is smaller during most of 
the evolution. However, due to the latter transition of w in Model B, 6e starts to decrease 
at a later time and its final value is larger then in Model A. 



3.3 Matter growth and ISW effect 

Once we know approximate analytical solutions to set the IC of the system of equations, we 
proceed to study the effects of EDE models on the linear evolution of matter perturbations 
and gravitational potential. 

In Fig. 5 we show the evolution of the normalised growth function D/a, where D = 
(a) /5m (1), and the logarithmic derivative of the growth function / = dlnD/dlna and its 
percentual difference against ACDM, 

A/ = 100 X - fACDM _ ^3_3g^ 

JKCDM 

In the ACDM model, both the normalised growth and / tend to a constant value for low 
a (high z) because D ^ a whenever DE is very subdominant. In this limit / — )• 1 and 
D/a — )• 1.34. In EDE models D/a and / vary for low a, but in the case of c^s = 0, DE 
perturbations tend to compensate the change in background evolution due to EDE and their 
values get closer to the ACDM values. In the case of Cgg = 1, DE perturbations are much 
smaller and practically do not compensate the change in the background evolution, thus the 
deviation from ACDM is more prominent. Besides this general behaviour of EDE models, 
we can see that Model B produces the most distinct evolution and variation between the 
two values of CcS- This is mainly caused by the rapid transition of its equation of state, 
which produces a very different background evolution and also allows DE perturbations to 
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Figure 5. Evolution of / (left top panel), its percentage difference against ACDM model (left bottom 
panel) and the normalised growth factor D/a (right panel), for a mode k = O.lhMpc"^. Line styles 
and colours for the different models are indicated in the legends. 



grow for a longer period when Cgfr = 0. Although the background evolution of Model B 
may be unrealistic, for instance, providing a rather low age of 13.01Gy, it is very valuable to 
clearly identify the changes that large DE perturbations may cause. Since current data on / 
have a precision around 10% [5], even in Model B, large DE perturbations would be weakly 
distinguished from the case of negligible DE fluctuations. However future surveys like Euclid, 
which may achieve about 1% precision [46], could make a much more significant distinction. 

The Sachs- Wolfe effect [47] causes the gravitational blueshift (redshift) of CMB photons 
when they fall in (escape from) potential wells near the time of last scattering. When the 
gravitational potential varies in time, the temperature variation of CMB photons is associ- 
ated with the Integrated Sachs- Wolfe effect (ISW). Assuming zero opacity the temperature 
fluctuation is proportional to: 



oc / (j)dri = (f) (r/doc) - 4> ivo) 



T 



(3.39) 



where r/dec and rjo are the conformal time at decoupling {zdec — 1100) and now, respectively. 
The ISW effect can be separated into two components: the early ISW, which is usually at- 
tributed to the influence of residual radiation energy after the redshift of last scattering and 
the late ISW, associated with the accelerated expansion at low redshift. Since we neglect 
radiation, we will only deal with ISW originated by DE. While ~ 1, the gravitational po- 
tential is nearly constant on all scales and the corresponding ISW during matter dominated 
era is negligible. However, for EDE models we have — 0.99 at high— z, so the background 
evolution alone should generate a distinct time evolution for (f>. Moreover, pressure pertur- 
bations in EDE directly source the time variation of the gravitational potential, Eq. (3.8), 
and it turns out that this contribution is more important. In order to visualise how EDE 
and its perturbations generate the ISW effect, in Fig. 6 we show the function 



F 



(3.40) 



so we have ^ oc </>(adec) (1 ~ -^(1))- 0i^6 can clearly see that the time variation of (p is 
greater in models with Ces = 1 than Ces = 0. Since the background evolution of Model A is 
similar to ACDM, in the case of Ccs = the ISW effect is also very similar to the ACDM 
model. Again, the general effect of models with null c^s is to compensate the changes of 
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Figure 6. Evolution of the function F, Eq. (3.40), as a function of the scale factor. The curve of 
Model A with Cgg = is almost indistinguishable from ACDM model. 

EDE in the background evolution. Although in principle ISW has limited power to constrain 
DE models due to its large cosmic variance, the observation of this effect is improving and a 
detection of 4.4(T level has been reported [48]. 

4 Nonlinear evolution 

Now we proceed to study the nonlinear evolution of DE and matter fluctuations, which can 
provide important information about the formation of DM halos. One difficulty that arises 
in this analysis is how to evolve DE fluctuations in the nonlinear regime. Many authors have 
addressed this issue using the Spherical Collapse Model, e.g., [42, 45, 49-51]. Independently 
of the approach used, the most important point we have to observe is when DE fluctuations 
actually become nonlinear, i.e., 5e ^ 1- The key property that controls this behaviour is 
the effective sound speed, c^s, which was shown can be constrained using future data on the 
abundance of galaxy clusters [43, 52]. 

As we showed in the linear regime, DE perturbations in models with non-negligible 
Ceff are at least few orders of magnitude smaller than 5^ on small scales. We also verified 
that, in this case, neglecting DE perturbations changes the value of the linearly evolved 6^ 
only about 0.02%. Moreover its value is proportional to the value of (j), which, before DE 
domination, is known to remain approximately constant even during the nonlinear regime 
[53, 54]. Hence we expect that the contribution of 6e for the gravitational potential will be 
very small during the nonlinear regime. For this reason, in the case of Ces = 1, we assume 
that DE perturbations can be neglected for the nonlinear evolution. 

For the case of Ces = 0, DE perturbations can have the same order of magnitude of 
matter perturbations, hence they can not be neglected, otherwise the linearly evolved (5m 
can change up to 20%. Evidently an intermediate behaviour should appear when Ccs > Cni, 
in this case one could incorporate the linear evolution DE perturbations with the nonlinear 
evolution of DM [42]. Here we want to focus on the two limiting cases: the one where 6e is 
most unimportant for the nonlinear evolution of DM, CgS = 1, and the other where de is the 
most important for it. 

We use the fluid approach to solve the SC model derived from the Pseudo-Newtonian 
Cosmology (PNC) [44, 51], which is suitable to treat the nonlinear evolution of perfect fluids 
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with relativistic pressure on small scales. The equations for a system with matter and DE 
fluctuations, with Ces = 0, can be written as: 



5m + e(l + <Jm) =0, (4.1) 

6e-3nwSe + e{l + w + 5e)=0, (4.2) 

a2 nnj2 

e + ne + - = — — {n^5^ + nj,) . (4.3) 



For the case of negligible DE perturbations we only need to evolve Eqs. (4.1) and (4.3), 
setting ile^e = in the latter, in this case we have the usual SC model for dark matter 
[55-57]. 

The most relevant quantity we compute using the SC model is the critical density 
contrast, 5c, defined by 

Sc {Zc) = SmL (zc) , (4.4) 

where 5mL is the linearly evolved matter density contrast with initial conditions such that 
non- linearly evolved 6^ has vertical asymptote at Zc, i.e., 

lim (Jjn (^) = oo . (4.5) 

2->2 + 

Numerically we need to determine some value above which we consider 5ra has diverged; we 
observed that choosing 10^ gives the classical value in Einstein-de-Sitter universe, Sc — 1.686, 
with an error of at most 0.5% for < z < 4. 

The function 6c depends both on the linear and nonlinear evolution. Although we do 
not have a framework to solve the relativistic SC model, we can compare the PNC and GR 
predictions for the linear evolution. It is clear that PNC will not give accurate results for 
scales of order of magnitude comparable with the horizon, however, for k = O.lhMpc"^, 
the values of 6^ given by PNC differ from GR at most 0.1% and decreases for smaller 
scales. Hence computing 6c using Eqs. (4.1) — (4.3) and its linearised version should be a 
good approximation. 

Now we need to provide the initial conditions to evolve the nonlinear equations. We 
choose values of 6m (oi) in order to get < Zc < 4. As we discussed in section 3.2, for the case 
of null Ceflf we can determine 6e (ai) using Eq. (3.30), use 6 (aj) = — (oi), assuming 6tn oc a. 
For the case of non-negligible Ces we assume 5m oc with p determined by Eq. (3.26). It 
is important to observe that all these relations can be derived from the PNC equations, 
Eqs. (4.1)-(4.3). 

We show the critical density contrast in function of z in Fig. 7. For the cases of Ceg = 
DE perturbations clearly make 6c closer to the ACDM values. This happens because, while 
at background level EDE changes the expansion rate and lowers Om, DE perturbations 
contribute to the gravitational potential, compensating the change in background. Note that 
the values of 6c that we found are consistent with small departures (less than 1%) from 
ACDM values [33, 58], even in the presence of large EDE fluctuations. 

The total mass of halos 

In models with clustering DE we must care about its contribution to the total mass of the 
halo. The presence of DE can either add or subtract mass from the forming DM halo. For the 
case of negligible Ceff , Ref. [45] argues that DE contribution is constant in time and evaluates 
it at the virialization radius, i?vir- 
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Figure 7. Evolution of Sc (zc) for the different models. Line styles and colour are shown in the legend 
box. 



The computation of i?vir depends on the details of the viriahzation process. In the 
usual SC model, in an EdS universe, the viriahzation radius is half of the turnaround radius, 
i?vir = -Rta/2. Oncc the viriahzation is reached one can define the DM overdensity of the 
formed halo. There are two common definitions: Av(-Zc) = {Pm{z^) + ^ Pmizv)) / Pm{zc) — 178, 
where is the redshift of viriahzation, and Av(-Zv) = (Pm(-2v) + 5 pm{z^)) / Pmizv) ~ 147, see 
Refs. [59, 60] for details and discussion. However when DE is present these values change in 
time and the viriahzation process also depends on the properties of DE [43, 45, 61, 62]. 

Here we will use the usual relation between turn-around and viriahzation radius -Rvir = 
i?ta/2. Although this is not valid in non-EdS models and the values of Ay can be considerably 
different, this simplified relation gives a good approximation for the fraction of DE mass to 
the DM mass: 

For instance, we checked that changing i?vir by 5%, which is common for inhomogeneous DE 
models Ref. [42], changes e by 1%. 

After solving the system of Eqs. (4.1) — (4.3) we find the value R using the conservation 
of DM mass: 

5„. + l = (5^, + l)(^§) . (4.7) 
The virialized mass associated with DM is then given by: 

Mm = 47r / ""^ dRR^{pm + 6pm) ■ (4.8) 
Jo 

Now we also have to define the DE mass contained in the same region where the DM halo has 
been formed, Mg. In Ref. [45] it is defined as the matter associated only with the fluctuations 
of DE, which we will call 

Mep = 47r / dRR^dpe . (4.9) 
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Figure 8. Evolution of e(z). Left panel: ACDM model and models A and B with Ccs = 1. Right panel: 
Model A only with perturbation contribution (dotted blue line) and total contribution (solid blue 
line), Model B only with perturbation contribution (dot-dashed red line) and with total contribution 
(dashed red line), all with Cgs — 0. 

However, this choice makes a discrimination between the DE and DM contributions. The 
latter is computed using its total density, Pm + Spm, whereas the former only considers the 
fluctuation contribution. Treating the two fluids on equal foot, the total contribution of DE 
mass can be evaluated based on the so-called active density in the Poisson equation in the 
presence of relativistic pressure, V^cp = 4-kG{p + 3p). Hence we have: 



The motivation to consider the total mass (background plus perturbative contributions) 
is that the astrophysical processes associated with galaxy clusters are sensitive to the total 
gravitational potential. This interpretation is also carried out by Ref. [63] in order to estimate 
the mass of the Local Group. Although this background contribution may be debatable we 
decided to evaluate it in order to verify whether it can produce a non- negligible effect. For 
models with Ces = 1, DE perturbations are much smaller than unity at the virialization time, 
hence they can be neglected in Eq. (4.10). 

In Fig. 8 we show the evolution of the quantity e for models with homogeneous DE (left 
panel), in which case we compute e using Eq. (4.10), and for models with inhomogeneous 
DE, computing e both with Eq. (4.9) and Eq. (4.10) (right panel). Note that due to the 
background contribution in Eq. (4.10), even in the case of Ces = 1 and ACDM models, DE 
can subtract only about 1% of DM mass. However this effect appears only at low redshifts, 
when ~ f^m- For the case of Ces = 0, the fraction of DE in a DM halo can be much larger, 
but now adding mass to the halo. In these models the DE contribution does not go to zero at 
higher z because of the combination of two factors: while fie decreases, but not to negligible 
values, w decreases as well, which in turn, according to Eq. (3.30), makes the ratio 6e/Sra 
grow. 

This is a characteristic behaviour of EDE models. As already shown in Ref. [45], for 
the perturbative contribution and constant w, around —1, e is at most 0.05 and tends to zero 
as z grows. In Ref. [43], for models with w = —0.8 and c^g = 10~^, the authors claim that 
maximal value of DE mass fraction in halos is 0.1%. However they assume DE fluctuations 
are small and can be treated within linear theory, which is not valid for the models with 
Ceff = we are studying. 




(4.10) 



5 Abundance of halos 



In this section we study the effect DE fluctuations on the abundance of halos. We decided 
to parametrize the mass function using the Sheth & Tormen (ST) prescription [64-66]: 



dn 
dM 




'M 

where a = 0.707, p = 0.3, A = 0.2162, 



1 r°° 

<^M = i^ j dkk^W^ {kR) P{k) (5.2) 

is the squared variance of the matter power spectrum, P(k), which we computed using 
the BBKS transfer function [67], smoothed with a top-hat window function, W{kR) = 
3{kR)~'^{sm{kR) — kRcos{kR)), where R is the scale enclosing the mass M = {An / 3) R'^ pm 
and is the comoving matter density. The ST mass function depends critically on the linear 
overdensity parameter 6c of dark matter, therefore also a small variation of this quantity will 
give a huge effect on the high-mass tail of the mass function. 

In order to clearly observe the effects of DE fluctuations we compute the number density 
of objects above a given mass at fixed redshift: 

We choose four different redshifts, namely z = 0,0.5,1,2. The actual number of halos also 
depend on a integral over the comoving volume, hence it depends both on background and 
perturbative properties of DE. For the models we consider the comoving volume is always 
smaller than in ACDM, see Fig. 2. 

We adopt as reference model the ACDM model with normalisation of the matter power 
spectrum ag = 0.776, in agreement with CMB measurements by the WMAP team [68, 69]. 
Since the background history and therefore the growth factor for the EDE models differ from 
the ACDM models, perturbations will evolve differently in the different classes of models. 
Therefore we decide to adopt the CMB normalisation: we fix the same amplitude of the 
perturbations at the CMB epoch and we rescale it by the ratio of the different growth 
factors. More quantitatively we have 

DACDuizdec) ^r- y,\ 

0'8,DE = 0-8,ACDM—f^ 7 ^ • i5-4j 

-t^DEl^decJ 

In Tab. 1 we show the normalisation of the power spectrum for the dark energy models 
considered in this work. 

As one can notice, the Model B is the one differing most from the ACDM model. 
The matter power spectrum normalisation when Cqs = 1 is 30% lower than the reference, 
therefore we expect that the mass function will be significantly different from the reference 
ACDM model. In Fig. 9 we show the ratio between the number density of dark matter 
halos for the EDE models and the ACDM model. We refer to the caption for the different 
line-styles and colours. 

It is immediately clear that for Model B structures are strongly suppressed at all red- 
shifts with respect to the ACDM model. Model B shows a lack of objects already at galactic 



-17- 



Table 1. Table of the normalisation of the matter power spectrum for the different models analysed. 
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Figure 9. Number density ratio between the EDE models and the reference ACDM model. The red 
dashed (blue short-dashed) curve shows results for Model A with CeS = (cog = 1), while the cyan 
dotted (orange dot-dashed) curve represents Model B for CeS = (ces = 1). The upper left (right) 
panel shows results for z = {z = 0.5), while the lower left (right) panel shows results for z = 1 
(z = 2). 



scales (M ~ 10^^ - 10^^ Mq/H) while at cluster scales (M w 10^^ Mq/H) the case with 
Ceff = 1 has basically no objects (see the upper left panel). At higher redshifts, the number 
of objects decreases considerably also at very low scales (see lower right panel for the results 
at z = 2). 

Model A suppresses structures less than Model B and at z = the Model A with Ceff = 1 
is quite similar to the Model B with Ceg = 0, differing by it at most of 10% for masses of 
the order of 10^^ Mq/H, while with the increase of the redshift, differences between the two 
tend to increase. The most similar model to the ACDM one is Model A with Ceff = 0. Also 
at very high masses differences are at most about 20%. It is also important to bear in mind 
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that the actual number of halos will also depend on the comoving volume, Fig. 2, which for 
EDE models is always smaller than in ACDM models. For low redshifts, z = and z = 0.5, 
the differences in the comoving volume is very small, but for z = 1 and z = 2 the actual 
number of halos in EDE models will be even smaller than the ratios in the lower panels of 
Fig. 9 suggest. 

It is also worth to note that models with Ces = have more objects than the case 
with Ceg = 1. This is easily explained in terms of the growth factor and of the evolution 
parameter 6c- Again the general effect of large DE fluctuations is to compensate the change 
in the background evolution by enhancing the gravitational attraction. 



The contribution of DE mass to the mass functions 

As already observed in Ref. [45], if DE can cluster it also contributes to the halo mass, thus 
we must compute the correction to the mass function due to this extra component. The 
actual contribution of DE crucially depends on whether it virializes and the time scale of 
this process. Moreover the merging history of halos formed at different times [70], which 
could contain different amounts of DE, should also be taken into account. A complete and 
accurate description of this corrections is a complex task, which depends on nature of DE, 
and is beyond the scope of this paper. Here we will compute a straightforward correction for 
a scenario where, in the case of Ces = 0, we assume that DE virializes together with DM on 
the same time scale. Hence, once the halo is formed, DE contribution is assumed to remain 
constant. This is probably the case in which DE fluctuations will mostly influence structure 
formation and consequently the mass function. 

Given a halo with DM mass M, its total mass will be M(l + e{z)), with e{z) defined in 
Eq. (4.6). Therefore, given the original DM mass function, j^^, we assume the corrected 
mass function is given by: 

We call the attention of the reader to the minus sign in M(l — e). Although the mass of 
a halo is changed by M — )• M(l + e) the use of this mass redefinition in the mass function 
would produce wrong results. For a positive e the halos become more massive than predicted 
by a model in which only DM clusters, hence more massive halos are expected. However, 
if one redefines the mass function using M(l + e), fewer massive halos would be predicted, 
which is just the opposite of what is expected. Therefore the natural correction of the mass 
function is the change of variable M — )• M(l — e). 

In Fig. 10 we show how much is the correction in the comoving density of objects above 
mass M relative to the same model without the correction, nc(> M)/n{> M). We present 
the results for models A and B, with Ceff = 0, only for the perturbative contribution (top 
panels). Although it is not clear whether the background energy density of DE is a stable 
contribution to the total mass of the halo, i.e., that is constant through the history of the 
halo, in order to have an idea of its influence on the abundance of halos, we also show the 
corrections due to the background contribution in the ACDM model (bottom panel). 

As can be seen in the these three cases, it is clear that even a small DE contribution for 
the halo mass can produce drastic changes at high mass. For Model B, the one with largest 
e, TLc can be many times larger then its version without the mass correction. In Model A 
the corrections are below 10% for small masses but can reach almost 60% for high masses 
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Figure 10. Ratio of comoving number density above mass M at fixed redshifts, computed with the 
mass function correction, Eq. (5.5), nc{> M) to the comoving number density in the same model 
without the correction n{> M). In Model B the corrections reach 60 times at z = 2. 



at high z. For ACDM, the corrections can be about 3% at z = and low masses, reaching 
more than 10% for very massive halos at high z. 

It is also clear that, although e gets smaller at higher z, the corrections do not necessarily 
diminish, as can be seen for the three cases with z = 2. Since the growth function, D, gets 
smaller for high z, the exponential decay of the mass functions is shifted to lower mass regions, 
hence even a small value of e at high z produces large modifications. The comprehension of 
this effect is very important because, despite the fact that very few massive galaxy clusters 
are expected at high z, the detection of a single massive distant galaxy cluster can be used 
to rule out DE models [71]. 

Now that we observed that the corrections in mass functions can be important, we have 
to ask ourselves whether they can change the behaviour that we previously found without 
taking into account such corrections, i.e., that all EDE models we consider have smaller 
density of objects then ACDM, Fig. 9. In Fig. 11 we plot the ratio of the corrected values 
of number density, ric, to the n in ACDM model. For small masses, the corrections indicate 
that our inhomogeneous EDE models actually present more objects than ACDM. However, 
nc(> M)/n{> M)/^cDM decreases with mass, except for Model A at z = 2. For Model B 
the suppression of structures due to its low as value is so strong that, even with the large 
corrections caused by its large values of e, at high masses it always has fewer objects density 
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Figure 11. Ratio of comoving number density above mass Af, computed witli the mass function 
correction Eq. (5.5), 7ic(> M) to the comoving number density in the ACDM model. 



than ACDM. 

In Model A a quite interesting behaviour can be observed. Since its and growth 
function are not very different than in ACDM, the correction due to DE mass turns out to 
be more important at high z and rapidly increases with mass. At 2: = 2 the number density 
in Model A is about 36% larger than in ACDM. Note that this effect occurs at high z and 
is strongly mass dependent. Ref. [71] points out that if unexpected massive clusters, within 
the ACDM or smooth quintessence paradigm, are present only at high z, DE clustering may 
not provide a consistent description because in these models cluster abundances are modified 
roughly by the same amount for both low and high redshifts. However, as we just have 
observed, the corrections on mass function due to DE mass may provide more abundant 
massive clusters at high z, without drastic changes for low z. Hence, if in the future the 
observation of a massive cluster falsifies ACDM and smooth quintessence models, it still can 
interpreted as an evidence of clustering DE. 

However, let us note once again that the actual number of objects also depends on the 
comoving volume. Therefore the behaviour observed in Fig. 11 can be modified, especially at 
high z. In order to verify what is the observable effect, we finally compute the total number 
of clusters above mass M and redshift z: 



N{> z,> M) 



dM'—— 
dM' 



M 



(5.6) 



In Fig. 12 we show the total number of halos above mass M and z = 1 (left panel) and 
z = 2 (right panel) for ACDM, Model A and Model B. The values for Model A and Model 
B are computed using the corrected mass function for the perturbative contribution only. 
For z > 1 the EDE models present fewer objects than ACDM. The same is true for z > 2, 
however Model A and ACDM are much more similar. Therefore, in these specific examples, 
we can see that the correction in the mass function due to DE mass is not enough to produce 
more massive and distant clusters than ACDM. Anyhow the main lesson we should take 
from these results is that if DE possesses fluctuations, its contribution for the halo mass can 
substantially change the abundance of massive galaxy clusters. The proper description of 
this effect is of major importance when testing inhomogeneous DE models with observations 
of such objects. 



- 21 - 




Figure 12. Total number of halos above mass AI and z = 1 (left panel) and z = 2 (right panel), 
Eq. (5.6). The solid black lines refers to ACDM model, the blue dashed (red dotted) lines to Model 
A (Model B), with Ccs ~ and correct mass function due to the perturbative contribution only. 

6 Conclusions 

In this paper we have studied the influence of inhomogeneous EDE both in Unear and nonlin- 
ear stages of structure formation. We have evaluated the matter growth, the ISW effect, the 
contribution of DE fluctuations for the total mass of the halos, the halo abundance relative 
to ACDM model and its corrections due the extra DE mass contribution. 

We have showed that the presence of large EDE fluctuations, i.e., models with Ces = 0, 
have the general property of making DM growth, ISW effect and halo abundance more similar 
to the predictions of the ACDM model than homogeneous EDE models, those with c^s = 1- 
In Model A, the one more similar to ACDM in the background evolution, the differences 
in / = din D/ din a against ACDM predictions are below 5% level, and the impact of DE 
fluctuations is even smaller. Eventually surveys like Euclid, which can provide data on / 
with 1% precision [46], will be able to distinguish between homogeneous and inhomogeneous 
EDE models. 

The analysis of the number density of halos initially showed that all EDE models provide 
fewer density of massive clusters, independently of the redshift considered. For models with 
nearly smooth EDE this conclusion remains valid. However, for inhomogeneous EDE models, 
once we account for the extra halo mass associated with DE fluctuations, which we showed 
that can be of the order of 10%, and make the corresponding correction in the mass function, 
this situation may change. We saw that, at high redshifts, the corrected number density can 
be larger than in ACDM. This was observed for our Model A with c^s = at z = 2. However, 
after computing the total number of halos, the impact of this correction is supplanted by the 
effect of the smaller comoving volume of EDE models and ACDM still presents more massive 
objects. 

It is important to stress that, in principle, the magnitude of DE fluctuations that we 
found and the corresponding impact on observables that we have studied is specific of EDE 
models. If DE has non- negligible energy density at intermediate and high redshifts the 
equation-of-state parameter w is close to 0, which, according with Eq. (3.30), enhances the 
magnitude of DE fluctuations. However we also showed that inhomogeneous EDE models 
actually make predictions more similar to ACDM than their homogeneous counterparts. 
Therefore we conclude that if the accelerated expansion is caused by an inhomogeneous EDE 
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model it will be challenging to distinguish it from the Cosmological Constant. 
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